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ABSTRACT We present an evolutionary model for "7-ray line emission from superbubbles 
based on shock acceleration of metal-rich stellar ejecta. Application to the Orion OBI association 
shows that 7-ray lines at the detection threshold of the SPI spectrometer aboard of INTEGRAL 
are expected, making this region an interesting target for studies of the interaction of supernova 
shocks with the interstellar medium. 
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1. INTRODUCTION 

The claim for a detection by COMPTEL (Bloemen et al. 1994) of an intense 
flux of 3-7 MeV gamma-rays from the Orion molecular complex, attributed most 
naturally to 12 C and 16 de-excitation lines, has led many authors to re-consider 
the nature and impact of energetic particles (EPs) in the interstellar medium (ISM). 
Although re-analysis of COMPTEL data suggests now that the observed emission 
was an instrumental artefact (Bloemen, these proceedings), the former "detection" 
raised the question about the possible existence of a low-energy, C and O enriched 
cosmic-ray component. Indeed, independent of the COMPTEL result, new ob- 
servations relating to the Be and B abundances in the early Galaxy support the 
existence of such a component (e.g. Gilmore et al. 1992; Duncan et al. 1992; Casse 
et al. 1995). A mechanism revived to explain low-energy C and O rich cosmic-rays 
has been the acceleration of particles in a superbubble resulting from the intense 
energetic activity of an OB association inside a molecular cloud (Bykov & Bloemen 
1994; Parizot 1998). Strong stellar winds and supernova (SN) explosions fill the 
superbubble with both energy and enriched material to be accelerated by the nu- 
merous secondary shocks and by magnetic turbulence resulting from the interaction 
of shock waves (from winds and SNe) with each other and with dense clumps in- 
side the bubble (Bykov & Toptygin 1990; Bykov & Fleishman 1992). The resulting 
energy spectrum is expected to be very hard (oc E^ 1 ) up to a cut-off energy, Eq, 
of ~ 100 MeV/n. As for the chemical composition of the EPs, it is clearly related 
to the composition of stellar winds and SN ejecta, although some contamination by 
swept-up and/or evaporated material is likely to occur. 

In this paper we calculate the 7-ray line emission associated with such a sce- 
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nario. As we believe that the Orion complex associated with the Orion-Eridanus 
supcrbubblc represents the most favoured target for a detection, we normalise our 
results to the distance of Orion (450 pc) and the stellar content as inferred from 
observations of the Orion OBI association (Brown et al. 1994). 

2. BASIC INGREDIENTS OF THE MODEL 

The first step in our model calculation consists in the evaluation of the enrich- 
ment of the superbubble by stellar winds and SN ejecta as a function of bubble age. 
For this purpose we follow the evolution of a coevally formed OB association, char- 
acterised by an IMF of slope T. The enrichment is calculated using the stellar yield 
compilation of Portinari et al. (1998) who combined the Padova stellar evolutionary 
models with SN models of Woosley & Weaver (1995). Additionally, yields for the 
production of radioactive 26 A1 have been taken from Meynet et al. (1997), Woosley 
& Weaver (1995), and Woosley et al. (1995). To determine the parameters of the 
superbubble "blown" by the association, we derive the time dependent mechanical 
luminosity of the OB association from the evolutionary tracks of the Padova group. 
Using this luminosity, we solve the dynamical equation for a spherical, homogeneous 
bubble (e.g. Shull & Saken 1995). The characteristic density and temperature of 
the bubble interior is dominated by the "mass loading" from evaporated gas off 
the shell. This mass loading dilutes the bubble interior with ambient ISM mate- 
rial which we assume to have solar composition. We calculate the conductive mass 
evaporation from the shell into the bubble by solving the equation of classical, un- 
saturated conductivity (e.g. Shull & Saken 1995). Even if we disposed of reliable 
stellar evolutionary tracks giving the composition of the winds and the SN ejecta, 
we would still have to evaluate the mixing of the ejecta with the evaporated ISM. 
To avoid such a hazardous attempt, we consider two extreme scenarios, in which 
the EPs are made of the stellar ejecta alone (models P), or a perfect mixture of the 
ejecta and the evaporated ambient material (models D). 

The second step of our calculation consists of accelerating the enriched material 
within the superbubble assuming a constant acceleration rate during a time tq 
following each SN explosion. The EP spectrum is thus normalised so that the 
energy injection rate E is equal to Esn/tq, where Esn = 10 51 erg is the SN energy. 
To calculate the time scale To, we assume that each new supernova influences and 
provides energy to a region of size L around its explosion site, in which particles are 
accelerated with an efficiency rj ~ 10~ 3 (Bykov and Fleishman 1992; Parizot 1998). 
Further assuming that the extension of the region in which particles are accelerated 
increases as L = v^t, where va — 200 km/s is the Alfven velocity, we find that the 
total energy injected in the form of EPs after time t = r is Eep = V n b^nv A TQ (E), 
where (E) is the mean EP energy, averaged over the assumed spectrum, and rib 
is the density of the superbubble interior. Equating -Eep to Esn, we obtain an 
estimate for t , which we then use to normalise the EP spectrum and thus the 
7-ray fluxes. For typical values of (E) = 100 MeV/n and rib = 10~ 2 cm~ 3 , we 
obtain r ~ 10 5 years, corresponding to an acceleration power of ~ 3 10 38 erg/s. As 
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FIGURE 1. Orion-Eridanus super bubble evolution assuming Y = —1.7 and 53 stars in the 
mass interval 4-15 Mq (Brown et al. 1994). The leading letter indicates if pure (P) or fully 
mixed (D) ejecta were accelerated. Subsequent letters indicate analytic (A) or Monte Carlo 
simulated (MC) stellar populations. 



argued by Bykov and Fleishman (1992), the energy spectrum of the EPs depends 
on their feedback over the magnetic turbulence and the shock waves system inside 
the bubble. Any detailed calculation of this spectrum being out of the scope of this 
paper, we consider here the cut-off energy, Bo, as a free parameter with values in 
the range 3 - 3000 MeV/n. 

3. APPLICATION TO ORION AND DISCUSSION 

The evolution of the Orion-Eridanus superbubble composition and the predicted 
7-ray line emission is summarised in Fig. 1. On the one hand we populated the IMF 
using Monte Carlo samples that are compatible with the present Orion population 
(Brown et al. 1994). On the other hand we studied the academic case of an 'analytic' 
stellar population where the IMF is densely populated by 'fractional' stars. While 
the latter case provides the average 7-ray line emission, the Monte Carlo sampling 
gives us the possible scatter around this average. 

Our models predict 4.44 and 6.13 MeV 7-ray line fluxes of the order of a few 
10~ 5 ph cm~ 2 s -1 , i.e. around the expected threshold of SPI for broadened lines 
(Jean 1996). We want to emphasise that this value is only an order of magnitude 
estimate due to the intrinsic uncertainties in our simplified model. Nevertheless, the 
result indicates that Orion is still an interesting target for the observation of 7-ray 
excitation lines due to its proximity and star formation activity. Our model predicts 
an 26 A1 production around 10 -4 Mq, corresponding to 1.809 MeV line fluxes of 
~ 610~ 6 phcm _2 s -1 . This is compatible with the upper limit of COMPTEL 
(Oberlack et al. 1995), and again is at the detection threshold of SPI. In particular, 
the observation of either the 1.809 MeV line or the excitation lines (or both) will 
severely constrain the model parameters and hence provide important information 
about shock induced particle acceleration. 

Among the most interesting observables is the 12 C*/ 16 0* line ratio. For ejecta 
mixed with the evaporated ISM, the ratio is always very close to the solar value 
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FIGURE 2. 7-ray line fluxes and 12 C*/ 16 0* line ratio as function of cut-off energy Eq for 
different association ages tg (= 10^ years; labelled by an additional '-tg'). 



(~ 0.76 for Eq = 100 MeV/n). For pure ejecta, this ratio may deviate significantly 
from the solar value, with values depending on the presence of a very massive star 
(M > 5OM ) in the association (like in simulation MCI). However, i2c*/ 16 0* is 
also very sensitive to the cut-off energy E due to the different energy dependencies 
of the excitation cross sections (cf. Fig. 2). Additionally, for Eq < 20 MeV/n the 
acceleration time scale To becomes too long, and hence the injection power too 
small, for significant 7-ray line emission. The ambiguity of interpreting a given 
line ratio from the 3D-space of parameters (Eq, association age, dilution) may be 
removed by jointly studying additional line ratios, e.g. LiBe*/ 16 0* (where LiBc* 
refers to the so-called Li-Be feature around 450 keV). We will discuss the expected 
correlations in a separate paper where we also give more detailed information about 
the modelling procedure (Parizot & Knodlseder 1999, in preparation). 
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